Testing gravity on kiloparsec scales with strong gravitational lenses 
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Modifications to GR generically predict time and scale-dependent effects which may be probed by 
observations of strong lensing by galaxies. Measurements of the stellar velocity dispersion determine 
the dynamical mass whereas measurements of the Einstein radius determine the lensing mass. In GR 
these two masses are equal; in alternative gravity theories they may not be. Using measurements of 
the stellar velocity dispersion and strong lensing around galaxies from the Sloan Lens ACS (SLACS) 
survey we place constraints on lensing in modified gravity theories and extend previous studies by 
applying this data to explore its dependence on various properties of the lens such as the lens redshift 
or mass and thereby constrain scalar-tensor, f(R) gravity theories, and generic parameterizations 
of deviations from GR. Besides applying the observations to these specific gravity theories, the data 
places a constraint on a generic dependence of modifications to GR on the lens mass and redshift. 
At the 68% confidence level we find that the ratio between the lensing and dynamical masses can 
only vary by less then 50% over a mass range for the lens galaxies of 10 12 < M/Mq < 10 14 and less 
than 40% over the redshift range 0.06 < z < 0.36. 

PACS numbers: 95.30.Sf, 04.50.Kd,04.80.Cc 



I. INTRODUCTION 

The ability to test our basic understanding of gravity 
has been surprisingly limited [l|. Most precision tests 
have concentrated on the motion of the planets and light 
within the solar system or the motion of binary pul- 
sars. Although measurements in the solar system have 
reached the level of testing deviations from general rela- 
tivity (GR) to one part in 10 5 [2j,|3|], they only constrain 
theories in the weak gravity limit, on scales of an AU 
(r ~ 10 12 cm), and at a single redshift, z = 0. Binary 
pulsar systems similarly test theories at z = and on rel- 
atively small scales but have the added aspect of testing 
gravity in the limit where it is large due to the compact 
nature of neutron stars [l|. 

The discovery of an accelerated cosmic expansion [3, Hi 
has led to a flurry of theoretical activity. Although it is 
possible to explain the accelerated expansion within GR 
by introducing a cosmological constant or new cosmic 
scalar field or other source of energy density, another 
approach is to consider these observations as the first 
observational indication of a need for modifications to 
Einstein's theory of GR [6( . 

Many groups have attempted to explain the accel- 
erated expansion within alternative theories of grav- 
ity. These observations can be explained as the result 
of the dynamics of a scalar field within a generalized 
scalar tensor theory 0413 ■ Another proposal modifies 
the Einstein-Hilbcrt action by the addition of a general 
function of the Ricci scalar, f(R) [ill . H3 |. Depend- 
ing on the functional form of the function f(R), such 
a term may give rise to late-time accelerated expansion. 
Refs. [TUEI, P ro P ose d a five-dimensional theory of grav- 
ity which may lead to an epoch of late time accelerated 



expansion. In addition to studies dedicated to the obser- 
vational consequences of specific modified gravity theo- 
ries there has also been interest in parameterizi ng g eneric 
deviations from GR on cosmological scales [l5T - [l8l ] . 

Many aspects of these modifications can be constrained 
or even ruled out by considering tests of gravity made in 
the solar system or through observations of binary pul- 
sars [19l422l |. However, given that these theories are nat- 
urally dynamical and scale dependent, solar system and 
pulsar tests can be of limited use. Therefore it is impor- 
tant to test gravity at a variety of scales and redshifts. 
In particular, since those modified gravity theories which 
are able to produce an epoch of late-time accelerated ex- 
pansion must become dynamically important when the 
acceleration starts to dominate (around z ~ 0.5 [231 ]) 
probes of modifications to GR around these redshifts are 
of the greatest interest. 

Observations of strong lensing around galaxies present 
an important and unique opportunity to probe modifica- 
tions to GR over a range of redshift and on kpc scales. 
The idea to use strong galaxy lenses to constrain modifi- 
cations to GR was first proposed in Ref. [24j . Following 
this, Ref. [25[ was the first to use a data set of 15 strong 
lenses from the Sloan Lens ACS survey (SLACS) [HH^. 
The main difference between the analysis presented here 
and the one in Ref. [25J is that here the analysis is ex- 
tended beyond constraining a universal value for 7ppn 
and uses the observations to place constraints on how 
7ppn may depend on various properties of the lens, such 
as its redshift or mass. This paper also extends the anal- 
ysis to the full SLACS data set of 53 lens systems as well 
as uses a more realistic model for the luminosity profile 
of the lens galaxy (a Hernquist profile, as opposed to a 
power- law profile used in Ref. [25|). The use of a more 
realistic luminosity profile leads to a significant shift in 



2 



the best fit 7ppn- 

This paper is organized as follows. In Sec.|TT]we present 
how scalar modifications to gravity affect the dynamics 
of massive test particles (i.e., stars) and photon trajec- 
tories differently. We discuss how a comparison between 
the dynamics and lensing signal leads to a test of gravity. 
In Sec. Mil we discuss the predictions from general scalar 
tensor theories with a massive scalar field. In Sec. II VI we 
discuss the predictions from f(R) gravity and empha- 
size its ability to rapidly suppress any deviations from 
GR and how this transition presents unique observational 
signatures depending on the mass of the lens galaxy. In 
Sec. IVII we discuss how measurements of strong galaxy 
lenses can be used to constrain modified gravity theories. 
We present our conclusions in Sec. IVIII 



a = j V ± (* + $)ci£, 
= 2p V±$ N d£, 



deflection of the image of a background source through 
an angle a is given by 



(5) 
(6) 



where Vj_ is the gradient transverse to the photon's un- 
perturbed trajectory and d£ is a length element along 
that trajectory and we have assumed that /j, is inde- 
pendent of position on the relevant scales. In terms of 
the bare potential, $tv, observations of stellar dynamics 
through the spherical Jeans equation [28[ measure the 
combination 



II. LENSING AND DYNAMICS IN 
WEAK-FIELD LIMIT OF MODIFIED GRAVITY 
THEORIES 

One of the basic ways that we can distinguish between 
different theories of gravity is through a comparison be- 
tween the predicted and observed motion of test particles. 
Such tests compare the motion of photons (which move 
on null geodesies) and the motion of non-relativistic mas- 
sive particles (which move on time-like geodesies). The 
different types of geodesies are sensitive to different com- 
ponents of the metric and hence their comparison allows 
us to measure those components. 

We start with the line element corresponding to weak 
gravity, 

ds 2 = -(1 + 2^)dt 2 + (1 - 2<$>)5 ij dx i dx j , (1) 

which has been written in the conformal Newtonian 
gauge and has introduced the two Newtonian potentials, 
^ and The 'bare' Newtonian potential is given by the 
usual Poisson equation, 



V 2 $A,(f) = AnGp. 



(2) 



A general modification to gravity introduces two new 
equations. One relates the potentials to the underlying 
mass density, 

ly 2 ($ + *) = iTTfiGp. (3) 
The other relates the potentials to one another 

-j- = 7ppn- (4) 

Note that GR is regained when p = 1 and 7ppn = 0. 
In general both p and 7ppn can depend on a variety of 
quantities that determine the space-time such as the local 
mass density, position, redshift, and so forth. We will see 
specific examples of this in the following sections. 

In order to distinguish between the two Newtonian po- 
tentials we compare the dynamics of stars within a galac- 
tic halo and the deflection of light around the halo. The 
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1 + 7ppn 



N ■ 



(7) 



Therefore with a knowledge of &n Eqs. © and {2} show 
that a comparison between lensing and stellar dynamics 
provides a measurement of 7ppn- 

In Appendix [A] we derive how a general scalar mod- 
ification to the GR field equations leads to a modified 
relationship between the lensing and dynamical masses. 
The effects of these modifications can be compactly writ- 
ten in terms of an effective source of stress energy that 
we denote T c g defined in Eq. (|A7|) . 



III. GRAVITATIONAL LENSING IN GENERAL 
SCALAR-TENSOR THEORIES 



Scalar-tensor theories of gravity [19|, I29hl32j present us 
with an example of a class of modified gravity theories 
that naturally predict a redshift-dependent 7ppn- 

We will consider a general scalar-tensor theory defined 
by the action 
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The gravitational field equation is given by 



(8) 



^bd(i) 
9 



]^{d a f) 2 + Uf + U{f) 



d^ipdvip - V^V u ip = kT^, (9) 



with the scalar field equation given by 



'■P 



(cM 2 +2f/', (10) 



where a prime denotes differentiation with respect to the 
field if. 

Linearizing the scalar field around its cosmological 
value, f — fo + 5f, and specializing to a static case 
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we have 



T cff = V 2 <V + 2U'6<p, 



W 2 6ip 



3 + 2cj B d(^o) 



m 2 Sip, 



(11) 
(12) 



where the scalar field mass, m, is a function of the back- 
ground field and derivatives of the potential U and Brans- 
Dicke function, wbd, whose exact form is not needed for 
this discussion. In general scalar-tensor models which 
produce late-time acceleration have both m ~ H and 
U ~ H so that the mass and potential are negligible on 
kpc scales. In the absence of a potential and scalar field 
mass we then find the standard result 13011 
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7ppn 



2uj B d[<Po(z)} ' 



(13) 



and 7ppn then depends on time through the background 
evolution of (p$. 

Soon after the expansion of the universe was shown 
to be accelerating many groups proposed scalar-tensor 
models as an explanation. In order to produce models 
with expansion histories in agreement with observations 
Refs. |7H10| established an algorithm by which observa- 
tions of both the expansion history as well as the growth 
of structure would enable a complete determination of 
the scalar-tensor theory. In particular, Refs. 0, H| real- 
ized that the specification of the expansion history [in 
the form of H(z)] allows a reconstruction of <p(z)\ in the 
absence of precise measurements of the growth rate a 
functional form for U(tp) must be specified from which 
follows 7ppn(z)- 

An interesting case considers the ability to produce 
accelerated expansion by introducing a cosmological con- 
stant, Cly, within a scalar-tensor theory which is less than 
its value required in the standard ACDM cosmology (i.e., 
Ha ~ 0.7). As discussed in Ref. [8| these models remain 
viable until some z max (which is typically of order unity) 
after which the theory becomes inconsistent 1 . There- 
fore, scalar-tensor theories with a cosmological constant 
fly < f^A can only explain the observed accelerated ex- 
pansion up to z m ax- This implies that an observation of 
the expansion history at redshifts greater than z max can 
rule these models out. Because of this Ref. Q emphasizes 
that, as opposed to solar system observations, measure- 
ments of the luminosity distance, Dl(z), at larger red- 
shifts will place the most stringent constraint on these 
theories. Here we shall see that a measurement of 7ppn 
at a range of redshifts less than z max can also serve to 
distinguish between these models and ACDM. 

The choice of f2y in scalar-tensor theories fully speci- 
fies the potential U which then allows a calculation of 
7ppn(^)- We show the evolution of 7ppn for various 



choices of fiy in Fig. [TJ Note that for the scalar-tensor 
models considered here the parameters have been chosen 
so as to be indistinguishable from GR at z — (i.e., in 
the solar system). Also note that as Qy approaches the 
ACDM value of 0.7 we regain GR at all redshifts and 
7ppn = 1- 



n v = o.6 
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FIG. 1: The evolution of 7ppn in scalar tensor theories which 
exactly mimic the flat ACDM expansion history with CIm = 
0.3. The various curves correspond to different choices of Qy 
going from zero on the bottom to 0.6 at the top in steps of 0.1. 
As Qy increases towards the value of Oa = 0.7 the differences 
between the scalar tensor theory and GR decreases. For all 
cases presented here the parameters of the theory have been 
chosen to be indistinguishable from general relativity at z = 
(i.e., in the solar system). 



IV. GRAVITATIONAL LENSING IN f(R) 
GRAVITY 

A sub-class of general scalar-tensor theories that 
present markedly unique predictions are those theories 
for which the Brans-Dicke parameter identically vanishes. 
Recently, these theories have been extensively studied in 
particular case of /(-R)-theories [111 [l2|. These theories 
contain a particularly interesting mechanism, known as 
the chameleon mechanism (33l . l34j , in which the modifi- 
cations to GR are rapidly suppressed around an object 
with sufficient density. This rapid change in the behavior 
of the theory presents a unique scale-dependent lensing 
signature. 

The action for /(i?)-theories takes the form 



1 As is discussed in detail in Ref. 
carries negative energy. 



for z > Zmax the graviton 



2k 



^[R + f{R)] + S„ 



(14) 
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where f(R) is a function of the Ricci scalar R and S m is 
the matter action. The gravitational field equation can 
be written as 



[1 + /'(£)] + [Rf'(R) - f + 2D/'(i?)] 



-V M V M /'(ii) = kT, 



(15) 



In this theory the Ricci scalar becomes a dynamical quan- 
tity whose equation of motion is determined by the trace 
of the field equation, 



□/'(#) = i(/vr _ /',./,')] -f 2/ 



(16) 



In the limit where / — >■ 0, Eq. (|16p implies the usual 
algebraic relationship R — —kT and GR is regained and 
T is the trace of the usual stress-energy tensor. 

Using the trace equation to rewrite the gravitational 
field equation this theory produces 



T, 



cff 



3 k 



[kT + R] 



3 k 



[2R f(R) - f] . (17) 



Solutions to the trace equation, Eq. (fToj) . determine the 
lensing predictions for this theory. To understand these 
solutions, we rewrite the trace of the field equation as 



Df'(R) 



dV 
df'(R) 



= 0, 



with 



(18) 



(19) 



We note that for functions f(R) which reproduce the ob- 
served expansion history, the minimum of this potential 
yields the general relativistic relationship between R and 

T, R = -kT m. 

Energetics drive the solution of Eq. (fT6"]) towards two 
limiting cases [35| . If f'(R) remains close to its asymp- 
totic, cosmological, value as we move within the galaxy 
this trades the energy cost of fixing the f'(R) at a high 
point in its effective potential against the gain in main- 
taining a nearly homogeneous field. In this case, the 
solution to Eq. (fT6|) can be found by linearizing f'(R) 
around its cosmological value and we have T e g = T/3 so 
that 7ppn = 1/2. On the other hand, if the scalar cur- 
vature is able to reach the minimum of its potential this 
will be at a cost in gradient energy since the scalar cur- 
vature will have to transition from its asymptotic value, 
Ro, to the general relativistic value R = — kT. At the 
minimum of the potential deviations from GR are highly 
suppressed. 

Within a given object, far away from the center the 
scalar curvature starts off near its asymptotic value, Ro, 
and evolves with radial distance from the center. If the 
object is too 'small', in a sense we will make clear in a 
moment, then R ~ Rq throughout the object and devia- 
tions from relativity will be of order unity. On the other 



hand, if the object is compact enough then the scalar 
curvature is forced to the minimum of its potential and 
R = —kT within some radius ro. Within that radius de- 
viations from relativity are highly suppressed and we say 
that the object is 'screened'. 
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FIG. 2: The combination K,p(r)r 2 for an NFW profile [Eq.ET] 
for masses incremented by an order of magnitude between 
10 10 Mq (bottom curve) to 10 13 M (top curve). The out- 
ermost point where |/'(i?o)| intersects these curves indicates 
the transition radius, ro, inside of which deviations from rela- 
tivity are suppressed. For a given value of \f'(Ro)\ those halos 
with a mass below some threshold will not be screened and 
will therefore exhibit order unity modifications to gravity. 

As shown in Refs. [22|, HH screening occurs within a 
radius ro implicitly given by 



\f'(Ro)\ < K P (r )r 2 , 



(20) 



where Rq is the value of the scalar curvature on cos- 
mological scales and p is the local value of the density. 
Note that if the density of an object is too small then 
Eq. ([20)1 is not satisfied at any radius and the object is 
completely unscreened with order unity deviations from 
GR throughout. One can think of |/'(i?o)| as determin- 
ing a characteristic gravitational potential for the theory 
so that when the local gravitational potential (npr 2 ) is 
larger than this characteristic value deviations from GR 
are highly suppressed. Numerical solutions show that the 
transition from R ~ Rq to Rq <C R = — up occurs over a 
relatively short length-scale (see, e.g., Fig. 10 in Ref. |35j) 
so we will approximate it by a step function. 

The most stringent observational constraint on f(R) 
gravity theories comes from the requirement that it pass 
solar system tests. Measurements of the motion of light 
in the solar system has placed the constraint 7ppn,o = 



1 + (2.1 ± 2.3) x 10" 5 HSi]. In order for f{R) gravity 
theories to pass solar system tests the theory must sup- 
press deviations from GR within our halo leading to the 
constraint |/'(i? )| < 10" 6 j35|. 

Measurements of lensing around other galaxies can also 
serve to constrain this theory. In particular, Eq. (|20|) 
shows that /(-R)-gravity predicts a lensing signal around 
galaxies which depends on halo mass. To see this, con- 
sider an NFW halo [37| of the form 



p(r) = p c S c — 



-1 / x -2 

r 

1 + — 



(21) 



where r s is the scale radius and p c is the critical density 
of the universe. The amplitude S c — (A/3)c 3 /[ln(l + c) — 
c/(l + c)] relates the concentration to the virial radius 
with an over density, A = 119. We also assume the mass- 
concentration relation 



9 



1 + z 



M 



.12 x 10 l2 h- l M, 



(22) 



where h is the Hubble parameter in units of 100 km/(s 
Mpc) [HIH]. This relation reduces the NFW profile to 
a one-parameter family which we take to be dependent 
on the virial mass, M. Since pnfw 7 " 2 r f° r r < r s 
and pnfw^ 2 oc t 2 f° r r > r s it is clear that the inner- 
most point at which deviations from GR are suppressed 
in /(i?)-gravity will occur at the scale radius r s . 

Looking at Eq. (|20p and the NFW density profile we 
can see that halos with masses which satisfy 



|/'CRo)| >2k Pc S c (M) 



(23) 



will not be screened 2 . Since S C (M) decreases with de- 
creasing M this sets an upper limit to the mass of halos 
which can be screened given a value for \f'(Ro)\- This de- 
pendence is shown as the solid line in Fig. [3] For masses 
below this threshold the theory deviates from GR by fac- 
tors of order unity and strong lensing around these halos 
have 7ppn = 1/2. Therefore, a measurement of 7ppn = 1 
in lower mass galaxies would place a more stringent con- 
straint on |/'(i?o)|- 

Finally, we must also take into account that \f'(Ro)\ 
depends on redshift. For f(R) models which produce 
late-time acceleration we have \f'(Ro)\ oc l/H", where 
n > so that \ f'(Ro)\ increases as the universe expands. 



For a galaxy with pr 



with m > (which will 



2 Strictly speaking this condition only applies to halos which are 
isolated. A strong lens galaxy which sits within a larger halo 
may be screened by the larger halo even though the mass of 
the lens halo is below the threshold given in Eq. I|23[ l 40]. The 
local environments of strong lens galaxies can be approximately 
determined using photometric data and shows that for the data 
considered in this paper (the SLACS survey) the richness of the 
lens systems is on average a few with very few close companions 
[41 1 indicating that most systems should be sufficiently isolated 
for our purposes here. 
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FIG. 3: The minimum mass for a halo to be screened given 
a value for |/'(i?o)|- Solar system tests require that the 
Milky Way's halo be screened. Since the halo has a mass 
~ 10 12 M e this implies that \f'{Ro)\ < 10" 6 Strong 
lensing measurements around galaxies with smaller masses 
will have 7ppn = 1/2. The dashed line and upper x- 
axis show how |/'(-Ro)| varies as a function of redshift with 
\f'(R )\(z = 0) = 10" 6 . In general, for f(R) models that 
produce late-time acceleration \f'(Ro)\ oc 1/H(z) n . In this 
figure we show the evolution of \ f'(Ro)\ for the model pre- 
sented in Ref. [3^|. The dotted red line indicates the way in 
which this plot should be read: at a given redshift (upper 
x-axis) a vertical line intersects a given f(R) model (dashed 
curve); a horizontal line then intersects the solid 'screening' 
line; from that point a vertical line drawn to the correspond- 
ing halo mass (lower x-axis) gives the mass above which the 
halo is screened at that redshift. 



be true in the outer regions of the galaxy to ensure that 
the galaxy has a finite mass) this causes the 7ppn = 1/2 
region to propagate inwards as time progresses so that 
/(-R)-gravity predicts a redshift dependent 7ppn as well. 
The dashed curves in Fig.[3]show the evolution of \f'(Ro)\ 
as a function of redshift for the specific f(R) model found 
in Ref. HI, 



f(R) 



c 1 (R/m 2 ) r 



V 



(24) 



c 2 (R/m 2 ) n - 

where m 2 ~ Hq and n is a free index and the ratio Ci/c 2 
is set by requiring the model have the same expansion 
history as in a ACDM model with Hm and fl\ 



c\ Mm 
— « 6— 



C-2 



0/ 



(25) 



Therefore these models have two free parameters which 
we choose to be n and \f'(Ro)\(z = 0). 
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V. STRONG LENSING PARAMETERIZED BY 
GRAVITATIONAL SLIP 

In order to test the predictions of GR we must com- 
pare predictions and observations with other theories of 
gravity. Besides looking at other theories on a case-by- 
case basis it is more useful to parameterize modifications 
to GR and constrain the value of those parameters. This 
approach has proven very sucessful when interpreting ob- 
servations of the effects of gravity in the solar system 
in the form of the parameterized post-Newtonian (PPN) 
formalism [l[ . In this formalism the relationship between 
different parts of the metric are parameterized by con- 
stant coefficients and observations of the motion of test 
particles (both massive and massless) measure the values 
of these coefficients. 

The PPN formalism relies on the presence of localized 
sources of stress-energy and so cannot be applied without 
change to a cosmological context. Some studies have at- 
tempted to articulate ways in which to extend the PPN 
formalism to cosmological observations • The ba- 

sic idea of all of the currently proposed parameterizations 
is that current cosmological observations using the cos- 
mic microwave background, weak lensing, and evolution 
of large-scale structure, are only sensitive to the scalar 
part of the metric [i.e., the two Newtonian potentials \I/ 
and $ in the metric given in Eq. ([T])] . Furthermore, since 
any modifications to GR that accounts for a phase of 
late-time accelerated expansion is significant for z < f 
the only perturbed fluid variables that are dynamically 
important are the matter density (5 m) and velocity per- 
turbations (9m)- Considering only those gravity theo- 
ries where stress-energy is conserved gives two evolution 
equations (the continuity and Euler equations). There- 
fore a generic modified gravity theory is defined when 
two gravitational field equations are specified. 

The two scalar potentials form a linear combination 
which is determined by a Poisson-like equation but with 
a time and space dependent Newton's 'constant' p(z,x) 



V 2 (A<I> + 5*) 
A + B 



4ir/j,(z, x)Gp, 



(26) 



and their ratio can be parameterized by another time and 
space dependent function given in Eq. Q and repeated 
here, 



7p PN (z,x) 



(27) 



In GR we have A = 1, B = 0, G(z,x) = G, and 
7ppn(Z) x) — 1- Particular modified gravity theories can 
then be parameterized by how the functions G(z,x) and 
7ppn(z, x) depend on scale and time [TBI. Il7|. 

One particular parameterization, first proposed in 
Ref. supposes j PPN (z,x) = [I + w' p DE /p M (z)}^ 1 = 
[1 + zuq/(1 + z) 3 ] -1 so that order unity modifications 
turn on around the transition from matter domination 
to dark energy domination. As described in more de- 
tail in Ref. [42j this parameterization further chooses to 



maintain the scalar part of the (0, i) component of the 
Einstein field equations leading to a time dependent p,. 
In the next Section we show how observations of strong 
lensing around galaxies over a range of redshifts are able 
to constrain the value of ujq. 



VI. 7P PN FROM MEASUREMENTS OF STRONG 
LENSES 

The original idea of measuring 7ppn from observations 



of strong lenses was first discussed in Ref. [24j and was 
first applied to data in Ref. (25[. A qualitative under- 
standing of how observations of strong lenses can yield a 
measurement of 7ppn can be understood by considering 
the following simplified example [24{ . In this example the 
lens density distribution is given by a singular isothermal 
sphere with the observed line-of-sight velocity dispersion 
<7 bs- The lens then produces an Einstein ring with a 
radius, Re, given by [43[ 



R E = 47TO-; 



obs 



f 



7ppn \ DlDls 
~2 J Dg 



(28) 



where Dx is the angular diameter distance to the lens 
(L), source (5), and between the lens and source (LS). 
The observed spectra give measurements of the source 
and lens redshifts as well as of the line-of-sight velocity 
dispersion [44]. The angular diameter distances are ob- 
tained by fixing a fiducial cosmology although the choice 
of cosmological parameters does not significantly impact 
the final result. The data then yields a measurement of 



7PPN 



2tt 



Ds Rec 2 
DlDls a 2 ohs 



(29) 



As we will now describe, in practice the problem is more 
complicated since the density profile of the lens cannot 
be described by such a simple model. 

As in the simplified model that was just discussed, in 
order to measure 7ppn we must relate the observed stellar 
velocity dispersion to the lensing observations. First note 
that in an analogy with Gauss' law the deflection angle 
in a modified gravity theory parameterized as in Eq. ([3]) 
and (U]) for a circularly symmetric lens depends on the 
enclosed mass as 



AD LS pGM(9) 

d 8 D L e 



(30) 



where is a unit vector projected on the sky centered 
at the lens, and M(9) is the projected mass enclosed 
within the angle 6. The lens equation j43j then relates the 
observed Einstein radius 6e = Re/Dl to the enclosed 
mass and p 



Ri = 4 



DlDls 
D s 



pGM(R E ). 



(31) 



It is useful to associate an effective 'lensing' velocity 
dispersion with each lens system. The effective velocity 



7 



dispersion, a\ c 
stein radius, 



is denned through the measured Ein- 



RE = 47T' 



^lens D L D LS 



D s 



(32) 



and the lens equation [Eq. ([3"T])] allows us to relate this 
to the projected mass within the Einstein radius 



M(R E ) = 4tt 2 ° lDlS 



We can therefore write 



c 2 fiG D s 



I^lGM(Re) 

7T 



Re 



(33) 



(34) 



This 'lensing' velocity dispersion, af ens , should not be 
confused with the observed stellar velocity dispersion, 
cr bs- In order to measure 7ppn, the equations of hydro- 
static equilibrium are used to relate the observed stellar 
velocity dispersion to M(Re). 

The analysis presented here extends the model due to 
Ref. [45| . The total mass density of the galaxy is modeled 
as a power law 



Pm(t) 



Pm,o [ — 



(35) 



The stellar component is well ht by a Hernquist profile 

El 



(36) 



where is the total stellar mass and r* is a scale radius 
which can be written in terms of the effective radius of 
an i? 1 / 4 luminosity profile as r* = i? c ff/1.8153. Finally, 
we allow for a non-zero anisotropy in the stellar velocity 
ellipsoid which is constant with radius, 



(37) 



The radial velocity dispersion is found by solving the 
spherically symmetric Jeans equation 28j for the stel- 
lar component. Projecting out the line of sight velocity 
dispersion and performing a weighted average over the 
luminosity profile within a circular aperture of projected 
radius Ra the observed stellar velocity dispersion is re- 
lated to the model parameters through the expression 



2 

°obs 



f 2fiGM E (R E ) (Re 



p-2 



* 1 w sec ) 1 



7r (1 + 7ppn)-Re V r„ 

(38) 

where a scc [47J is the seeing (i.e., blurring due to atmo- 
spheric distortions), and g(p,/3, R^fr*, <r S cc) can be writ- 
ten in terms of integrals over hypergeometric functions. 
This 
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FIG. 4: The degeneracy between 7ppn and the slope of the 
total matter density p. In order to generate this curve we 
have used the mean values from the SLACS survey: Re = 
R e sf2 = 5 kpc with z L = 0.1 and /3 = . 



The ratio between the observed stellar velocity disper- 



sion and the lens velocity dispersion, f% = o-^ bs /a 
given by 



2 

lens 



7PPN 



Re 



p-2 



g(j>,/3,RA/r*,a s , 



(39) 



Observations of the stellar velocity dispersion and the 
Einstein radius yield measurements of f a . However, from 
Fig. @] it is clear there is a degeneracy between the slope 
of the density profile, p, and 7ppn- In order to make 
progress in measuring 7ppn a prior must be placed on 
p. It is important that any approach taken to place this 
prior be independent of the theory of gravity. One ap- 
proach, discussed in Ref. [Il], uses an assumed scaling 
law (related to the fundamental plane [49j]) that relates 
the power-law slope, p, to the luminosity, Einstein radius, 
and effective radius is used to estimate the average value 
for p within a population of lenses. In this approach a 
universal value for 7ppn would appear as a constant offset 
and does not affect the estimate of (p) = 1.959 ± 0.077 
[48| . However, we are interested in exploring whether 
7ppn may be non-universal. Furthermore, the analysis 
in Ref. [48| is done in the limit where the dark matter 
fraction of in the system is 1. This is a significant sim- 
plification since the average dark matter fraction within 
Re of a subset of these systems has been estimated to be 
25% 

Another approach is to place a prior on p given by its 
distribution measured in low redshift early-type galaxies 
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FIG. 5: The mean (left) and dispersion (right) for the distribution of v2/ct/\/1 + 7ppn as a function of the observed effective 
radius, i? e ff, and the observed Einstein radius, Re- The intrinsic distribution for the slope of the total density profile, p, and 
the velocity anisotropy, /3, were assumed to be equal to their intrinsic distributions measured in early-type galaxies at low 
redshift where more detailed measurements of the stellar kinematics can be measured. 



where more detailed kinematic data can be used to deter- 
mine the full density profile. This is the approach taken 
here and gives Op) = 1.93, a p = 0.08 and (0) = 0.18 
and = 0.13 [U, [HI- But as is clear from Fig. [U the 
constraint on 7ppn is very sensitive to the assumed (p) 
(see Fig. @|. To remove this dependence on (p) in our 
final constraints we marginalize over the average of value 
7ppn within the sample. Therefore, our final constraints 
come from the lack of any significant correlation between 
7ppn and properties of the lens such as its redshift and 
mass. 

Using Eq. (131))) the resulting distribution for 
V^/ct/VI + 7ppn was calculated on a grid of val- 
ues for 0.2" < R cff < 7" and 0.2" < R E < 4" using 
an aperture radius Ra = 1.5" and seeing a scc = 0.64" 
50] and is shown in Fig. [5] For a value of 7ppn, an 
underlying probability distribution for p and f3, and 
measured values for Re and i? e ff, Eq. ([39]) gives the 
underlying probability distribution of / CT , denoted by 
(dP I dfa) (7ppn) • The likelihood for 7ppn is then given 
by 



£( 7 : 



PPNj = 



(40) 



dP 



~Tf~ (7ppn) [U\Re, R eS }G(f a ;a u , f a )df a 



where G(f a \ <jf a , /°) is a Gaussian with a mean equal to 
the measured value of /° and a standard deviation, cr* , 
equal to the observational error on f a . 

Since the resulting probability distribution of 
V^/cr/Vl + 7ppn for fixed R e g and Re is well approxi- 
mated by a Gaussian it is straight-forward to derive the 
distribution of f a for any value of 7ppn, -Rcff, and Re- 



We first consider the case where 7ppn is a universal 
constant and the results are shown in Fig. [5] In that fig- 
ure the solid black curve shows the likelihood for the full 
SLACS survey (53 systems) The full joint likelihood 
gives the result 7ppn = 0.88 ±0.05 at 68% c.l. As we dis- 
cuss in more detail below, even though this result seems 
to be in conflict with GR, it is more of a reflection of 
a difference between the mean slope of the mass density 
of the galaxies in the local universe as compared to the 
SLACS galaxies. The other two dashed-dot curves are 
the likelihoods when considering a subset of the full sur- 
vey. The dot-dashed black curve on the right is the like- 
lihood for the original 15 systems considered in Ref. [25[ 
and the dot-dashed red curve on the left uses 14 of those 
systems but reflects differences in the mass modeling, sur- 
face brightness measurements, and measurements of the 
velocity dispersion. As discussed in Ref. [53| the differ- 
ence in mass estimates, measured surface brightness and 
velocity dispersion, provides a more realistic sense of the 
measurement errors for these quantities. The dominant 
effect on the measurement of 7ppn is the error in the 
measured stellar velocity dispersion. The analysis pre- 



sented here was slightly more conservative than Ref. 53 1 
and took the statistical errors in as quoted in Ref. 53 1 
while enforcing a minimum of 7%. The agreement, at the 
lcr level, between the original analysis (dot-dashed black) 
and the modified analysis (dot-dashed red) indicates the 
error budget used in this work appropriately incorporates 
errors in the mass modeling, surface brightness measure- 
ments, and velocity dispersion. 

Comparing the constraints presented here to the re- 
sults found in Ref. [25], 7 PPN = 0.98 ± 0.07, we find 
a similar error but a significantly lower best fit value. 
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FIG. 6: The joint likelihood for a universal 7ppn- The solid 
black curve is the joint likelihood for the full SLACS sample 
(53 systems) [k| which gives 7ppn = 0.88 ± 0.05; the black 
dashed-dot curve on the right is uses the original 15 systems 
considered in Ref. [25|] which gives 7ppn = 0.93 ±0.1; the red 
dashed-dot curve on the left is a re-analysis of the 14 systems 
in common between the original 15 and the full 53 systems 
which gives 7ppn = 0.81 ± 0.1. As described in Ref. [53J the 
re-analysis of these 14 common systems reflects differences 
in the mass modeling, surface brightness measurements, and 
measurements of the velocity dispersion. The effect these dif- 
ferences have on the inferred constraint to 7ppn is well within 
la and so is well described by the errors we have included in 
the calculation of the likelihood (which are dominated by er- 
rors on the measured velocity dispersion). All quoted errors 
are at 68% c.l. 



This difference can be explained by noting that an anal- 
ysis of the SLACS lenses assuming GR (i.e., 7ppn = 1) 
found an intrinsic distribution for the slope of the mat- 
ter density, p, of (p) = 2.0 and a v = 0.12 [HHj. Since 
the mean of this distribution is approximately 0.1 away 
from the mean of the low-z distribution used to constrain 
7ppn (i-e., (p) = 1.93) and given that the degeneracy be- 
tween p and 7ppn is nearly linear (see Fig. 2]) it follows 
that the best- fit value should be about 0.1 away from 
7ppn = 1- The main difference between the analysis pre- 
sented here and the one in Ref. [25| is that this analysis 
uses a more realistic model for the luminosity profile of 
the lens galaxy (a Hernquist profile), which was also used 
to measure p in the original SLACS data [5Cj , as opposed 
to Ref. [25| which used a power-law profile. 

The degeneracy shown in Fig. @] has a significant ef- 
fect when applying the low redshift measurements of the 
slope of the total mass density p to the SLACS objects 
and leads to an apparent disagreement with GR at the 
2.4cr level and it is clear that this constraint on 7ppn is 




0.1 0.2 0.3 

Z L 

FIG. 7: The measurement of 7ppn in individual lens systems 
as a function of the virial mass of the lens halo (top) and lens 
redshift (bottom). 

of limited interest. In particular, the lack of knowledge 
of the mean value of the slope of the total mass density 
in the lens galaxies leads to a bias in any constraint on a 
universal value of 7ppn- In order to remove this bias we 
marginalize over the mean value of 7ppn for the entire 
sample of lenses. In this way the data does not constrain 
a universal (constant) value of 7ppn but instead gives ro- 
bust constraints on how 7ppn may depend on properties 
of the lens such as its redshift and mass. As we saw in 
previous sections, various theories of gravity give specific 
predictions on how 7ppn depends on both lens redshift 
and mass. 

In order to demonstrate how a variable 7ppn is con- 
strained we applied the SLACS data to the two modified 
gravity theories described in Sec. IHII and IIV1 A variable 
7ppn is constrained by calculating the likelihood 

£oce-* 2/2 , (41) 

with 

X 2 = E { - A lPPN (z Ll , Mi)] /<7 7PPNii } 2 , (42) 

i 

where 7pp s N ; is the observed value, ct 7ppn , is the error 
on the observation, 7ppn(-2l,i, Mi) is the predicted value 
which depends on either the lens redshift or its mass, and 
A is an amplitude which takes values between 0.8 and 
1.2. Marginalizing over A (which is similar to marginaliz- 
ing over the bias when using measurements of large-scale 
galaxy clustering) allows us remove the degeneracy be- 
tween (p) and 7ppn (which, as discussed before, is nearly 
linear; see Fig. [4]) and to explore constraints which arise 
due to (the lack of) any correlation or observed relation- 
ship between 7ppn and lens redshift or halo mass. 
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The lens redshift is a directly observable quantity 
whereas the mass is not. To estimate the total virial 
mass, M, from observations we use an approximate rela- 
tion with the observed half-light radius R e g- Using the 
fact that the stellar to virial mass ratio is approximately 
0.01 and given the relation between the stellar mass and 
effective radius found in Ref. [54[ for elliptical galaxies 
we can write 



M = 7.5x 10 11 (J^j 



1.78 



Mr. 



(43) 



Fig. [7] shows the measured value of 7ppn as a function 
of both halo mass (top panel) and lens redshift (bottom 
panel). Note that the mean estimated virial mass is 11 x 
10 12 Mq which compares well with the mean virial mass 
determined through weak lensing (M ) = 14+g x 10 12 Mq 
which was made using a subset of the full 53 systems 
considered here [Hj]. 

We may approximate a generic non-universal 7ppn as 
depending linearly on some property of the lens. In par- 
ticular, letting x denote a property of the lens (such as 
its redshift) we model this generic dependence as 



7ppn(z) = 7o (mx) + m. 



x 



(44) 



where 7o(to, c ) is a constant which depends on the slope 
m x and Ax is the range of x over which we have observa- 
tions (i.e., the range of lens redshifts for a given survey). 
In order to remove any sensitivity to the mean value of 
7ppn we define 



7o(m x ) = 7ppn - 



m x x 
Ax' 



(45) 



where 7ppn is the mean value of 7ppn for all of the lenses 
in the survey. Here we will only be interested on con- 
straining how 7ppn may depend on the lens mass and 
redshift. In these cases, at the 68% confidence level (c.L), 
we find that m z = 0.15 ± 0.24 (0.06 < z < 0.36) and 
m M = 0.13 ± 0.36 (10 12 < M/M Q < 10 14 ). 



A. Constraints to scalar tensor gravity 

As described in Sec. IIIII scalar-tensor gravity generi- 
cally predicts a redshift dependent 7ppn- For the partic- 
ular case where we account for the observed accelerated 
expansion using a scalar-tensor theory with a cosmolog- 
ical constant fly < the gravitational lensing obser- 
vations place constraints on the value of Qy as shown in 
the left panel of Fig. [8] As Sly is made smaller the value 
of 7ppn varies more with redshift leading to tension with 
the SLACS data (see Fig. [T]) The data places the con- 
straint fl v > 0.3 at 68% c.L and Qy > 0.06 at 95% c.L 
Note that since we have marginalized over the ensemble 
average for 7ppn these constraints rely solely on the lack 
of any significant correlation between 7ppn and the lens 
redshift. 



Since these models are constructed to have negligible 
deviations from GR today measurements of either the 
expansion history or lensing and dynamics at high z pro- 
vide the only data which can place meaningful constrains 
on these theories. 



B. Constraints to f(R) gravity 

As described in Sec. IIVI f(R) gravity generically pre- 
dicts a lensing signal which depends on the mass of lens 
halo. Using the relationship between the mass thresh- 
old (below which 7ppn = 1/2 and above which is equal 
to unity) shown in Fig. [3] the SLACS data places a con- 
straint on the allowed values of \f'(Ro)\(z — 0). Con- 
straints to the particular model described in Eq. ([24]) 
are shown in the center panel in Fig. [H Placing the 
prior \f'(R )\(z = 0) < 10~ 5 for the particular case 
where n = 1 we find \f'(R )\(z = 0) < 1.8 x 10~ 6 
68% (2.5 x 10~ 6 95%); n = 3 yields \f(R )\(z = 
0) < 2 x 10~ 6 68% (2.8 x 10~ 6 95%); n = 5 yields 
|/'0Ro)IO = 0) < 2.2 x 10~ 6 68% (3 x 10~ 6 95%). The 
constraint to \ J'(Rq)\(z = 0) is weakened as n gets larger 
since, as demonstrated in Fig. [3j the larger n is the more 
\f'(Ro)\(z) decreases with increasing redshift. Therefore, 
for the same value of \f'(R )\(z = 0) a model with a 
larger n will have a smaller value of \f'(Ro)\(z ~ 0.2) 
leading to a weaker constraint at the SLACS redshifts. 

The constraints to f(R) gravity from the SLACS data 
are about a factor of 10 worse than constraints from 
solar system tests. However, they are several orders 
of magnitude better than constraints using observations 
of anisotropics in the CMB and measurements of the 
matter power spectrum (H^ . A survey which could 
measure lensing and dynamics around low-mass galax- 
ies (M < 10 12 Mq) could potentially place constraints 
on f(R) gravity which would improve upon solar system 
tests. 

Note that since we have marginalized over the mean 
value of 7ppn these constraints rely solely on the lack 
of any significant correlation between 7ppn and the lens 
mass. 

C. Constraints to the gravitational slip 

As described in Sec. |V]a particular way to parameter- 
ize deviations from GR is to introduce two new time and 
space dependent functions: a modified Newton's 'con- 
stant' fj,(z,x) and jppa(z,ot). Assuming that the modifi- 
cations of GR become important at the same time as the 
expansion starts to accelerate inspires the parameteriza- 
tion [H 



7ppn(z) 



1 



1+TU (1 + Z)- 3 ' 



(46) 



and measurements of the SLACS lenses places a con- 
straint on vjq as seen in the right-most panel in 
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FIG. 8: Constraints to modified gravity theories and parameterizations using the SLACS data. Left: The left panel shows 
constraints to the energy density in a cosmological constant in units of the critical energy density, Qy, in scalar tensor theories. 
As discussed in Sec. 1 1 1 1 1 scalar tensor gravity can explain the observed accelerated expansion with a value for the cosmological 
constant which is less than the value required in GR JIa- As Qy is made smaller the value of 7ppn varies more with redshift 
leading to tension with the SLACS data (see Fig. [TJ. The SLACS measurements place the constraint Qy > 0.3 at 68% c.l. 
and £ly > 0.06 at 95% c.l. Center: The center panel shows constraints to the value of |/'(f?o)| using the model presented in 
Eq. IHl). The three lines correspond to n = 1 (solid red) \f'(R )\{z = 0) < 1.8 x 10" 6 68% (2.5 x 10" 6 95%); n = 3 (dashed 
black) \f'(Ro)\(z = 0) < 2 x 10~ 6 68% (2.8 x 10~ 6 95%); n = 5 (dot-dahsed blue) \f'(R )\(z = 0) < 2.2 x 10 -6 68% (3 x 10" 6 
95%). Right: The right panel shows constraints to the 'gravitational slip' zuo which parameterizes the evolution of 7ppn in 
time as discussed in Sec.[V| The SLACS data places the constraints vjq = 0.25lg 27 lo 48 (68%, 95% c.l.) which is as restrictive 
as constraints derived from the cosmic microwave background, weak lensing, and evolution of large-scale structure [42| . 



Fig. [5] The SLACS data places the constraints -hjq = 

0.25toir -oil ( 68% > 95% C - L ) which is as restrictive 
as constraints derived from the cosmic microwave back- 
ground, weak lensing, and evolution of large-scale struc- 
ture 

Note that, as in the case of f(R) gravity, since we 
have marginalized over the mean value of 7ppn these 
constraints rely solely on the lack of any significant cor- 
relation between 7ppn and the lens redshift. 



VII. CONCLUSIONS 

Constraints to modifications of GR from data taken 
within the solar system or from binary neutron star sys- 
tems within our galaxy are very precise. For instance, 
radar ranging to the Cassini spacecraft leads to a con- 
straint 7ppn,o = 1 + (2.1 ± 2.3) x 1CT 5 [2]. Although 
these measurements place important constraints recent 
interest in modified gravity theories which can account 
for the observed accelerated expansion have focused at- 
tention on models with modifications that evolve with 
both time and environment. 

In the case of scalar-tensor theories of gravity time 
evolution is a natural consequence of introducing a new 
scalar degree of freedom. In the case of f(R) gravity the 
chameleon mechanism, where modifications to GR are 
suppressed in regions of high mass density, is a natural 
consequence of the fourth order nature of the modified 
field equations. Other theories, such as DGP gravity [HI, 
03] and the recently proposed Galileon [57| , also predict a 



non-universal 7ppn- It is therefore important to not only 
look for ways to constrain 7ppn in the local universe but 
to also investigate whether it may change depending on 
time, scale, mass, local environment, and so forth. 

Using measurements of stellar velocity dispersions and 
strong lensing around early-type galaxies from the full 
SLACS survey we have presented constraints to alterna- 
tive gravity theories which can account for the observed 
late-time acceleration. This analysis updates the results 
presented in Ref. [25| by including more realistic model- 
ing of the stellar component as well as by using the full 
53 systems in the SLACS survey. We also extended the 
analysis beyond constraining a universal value for 7ppn 
and applied the data to constrain 7ppn's dependence on 
the mass and redshift of the lens. 

Constraints to a universal value for 7ppn must be used 
with caution given the significant degeneracy between the 
slope of the total matter density, p, and 7ppn- The two 
methods discussed in the text which estimate (p) inde- 
pendently of the theory of gravity both have unquali- 
fied systematic errors. In the case where a scaling law 
is assumed the scaling law itself may introduce biases; 
the application of low-redshift observations to the higher 
redshift SLACS lenses may not be appropriate given that 
the structure of the galaxies may significantly evolve with 
redshift. To remove this uncertainty the analysis pre- 
sented here marginalized over the mean of the sample 
leading to constraints which only depend on how 7ppn 
correlates with lens redshift and mass. 

Attempts to constrain modifications to GR on Mpc 
scales using observations of galaxy clusters, weak lens- 
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ing, and galaxy surveys are complementary to the pre- 
sented here j5a - [62j . Observations of galaxy clusters allow 
a measurement of 7ppn though a comparison between 
the X-ray temperature or virial mass and measurements 
of strong lensing. However, these observations have a 
limited statistical significance leading to constraints on 
7ppn to ~ 50% [63[ ■ This is partly due to the lack of 
large homogeneous samples of clusters. It is also related 
to the fact that cluster dynamics are harder to model 
leading to larger systematic errors. 

Although the SLACS survey presents us with an op- 
portunity to constrain the non-universality of 7ppn the 
dynamical range of the the SLACS survey are limited: 
their redshifts range 0.05 < z < 0.35 and their masses are 
of order 10 13 M Q . Future surveys may be able to extend 
this to a higher redshift as well as to lower mass galax- 
ies [flij which would improve the results presented here. 
For instance, measurements of lensing around galaxies 
with masses M < 10 12 M Q could potentially place a 
more stringent constraint on f(R) gravity theories than 
solar system tests. Of course the challenge to measuring 
strong lensing around less massive galaxies is that the 
lensing cross section decreases with decreasing mass. 



Acknowledgments 

After completing this paper, we became aware of simi- 
lar work by Schwab, Bolton, and Rappaport. The author 
thanks Adam Bolton, Robert Caldwell, Daniel Grin, and 
Eric Linder for a careful reading of a previous version of 
the manuscript, and Kevin Bundy, Bhuvnesh Jain and 
Marc Kamionkowski for useful conversations. This work 
was supported by the Berkeley Center for Cosmological 
Physics. The author gratefully acknowledges the hospi- 
tality of the Aspen Center for Physics where some of this 
work was completed. 



Appendix A: Lensing and dynamics in weak-field 
limit of modified gravity theories 

We start with the line element corresponding to weak 
gravity, 

ds 2 = -(1 + 2V)dt 2 + (1 - 2<f>)5 ij dx i dx j . (Al) 
Einstein gravity is determined through the field equations 

(A2) 

where, writing the metric as g^ v = rjp V 



G py — ^FfiW) 



hp U we have 



G --- 



_ _ d M cr„ + rj^dpO-f , (A3) 



where 7 MI/ = hp V ~ (l/2)i]^h, h = h^r}^ , and 07, = 
<9 r 7 M r. Considering a scalar modification of the field 
equations to linear order (appropriate for several mod- 
ified gravity theories that lead to a late-time accelerated 



expansion 04l3, we wri te Gp U (i~ 1 Gp U + r]p U ai + 
a 2dfj,dv a 3 + d^az,d v }a.Q, where a 2 and 04 are evaluated 
on the background which may be time-dependent. Now 
we specialize to a gauge in which we set 

<Jp + dpF = 0, (A4) 

with F satisfying the equation 

dpd v F = a 2 dpd u a3 + d^a 5 d^ae- (A5) 

With this gauge condition we can put the linearized field 
equation in the form 

□7/w = _2k M [Tp.v - -^Tesn^uj , (A6) 



where we have defined 



2ai + a^Da^ 



(A7) 



and we can neglect d p a^d p ae since it will be second order 
in the perturbation. We take the stress energy tensor to 
be dominated by pressure-less matter so that 



(A8) 



Specializing to a static source the solution to the field 
equation which is asymptotically flat becomes 

7/** = -4/x {$jvW"0 + , (A9) 

with the 'bare' potential given by 



$n(x) = - — 



" 1 " { -' : ' ] dV 



and 



8ir J \x — x' 



<M2) = / J^rdV 



8tt J \x — x'\ 



(A10) 



(All) 



This new scalar degree of freedom will be determined 
through a field equation whose solution will depend on 
boundary conditions. If we suppose that the modifica- 
tion introduces a new mass scale, m, then the effective 
potential may depend on $ c ff ($iv, mr, z, VIm, ■ •■)■ 
Taking the line-element written in Eq. (1) we have 



so that 



We then have 



/i = 2/i(*-3$), 
7oo = -/i(3$ + *), 



* = H ($jv + $eff) 



(A12) 



(A13) 
(A14) 



(A15) 
(A16) 
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